Induced Current and Aharonov-Bohm Effect in Graphene 
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The effect of vacuum polarization in the field of an infinitesimally thin solenoid at distances much 
larger than the radius of solenoid is investigated. The induced charge density and induced current 
are calculated. Though the induced charge density turned out to be zero, the induced current is 
finite periodical function of the magnetic flux $. The expression for this function is found exactly in 
a value of the flux. The induced current is equal to zero at the integer values of $/$o as well as at 
half- integer values of this ratio, where $o = 2-nhc/e is the elementary magnetic flux. The latter is a 
consequence of the Furry theorem and periodicity of the induced current with respect to magnetic 
flux. As an example we consider the graphene in the field of solenoid perpendicular to the plane of 
a sample. 

PACS numbers: 81.05.Uw, 73.43.Cd 
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The Aharonov-Bohm effect [l[ , scattering of a charged 
particle off an infinitesimally thin solenoid, which is ab- 
sent in classical electrodynamics, has been investigated in 
numerous papers, see Review [2|. Both non-relativistic 
[l| and relativistic [!, 0, 0, Q equations have been consid- 
ered. Similar effects having topological origin have been 
studied in quantum field theory in 0, @]. Intensive in- 
vestigation of the topological effects in condensed-matter 
systems has been performed recently both experimen- 
tally and theoretically in 0, [l(J 11 , 12| • New possibilities 
to study topological effects in Quantum Electrodynam- 
ics (QED) have appeared after recent successful fabrica- 
tion of a monolayer graphite (graphene), see Ref. [IH 
and recent Review [14j. The single electron dynamics in 
graphene is described by a massless two-component Dirac 
equation 15, 17, [lj[ so that graphene represents a 



peculiar two-dimensional (2D) version of massless QED. 
This version is essentially simpler than conventional QED 
because effects of retardation are absent in graphene. 
However, the "fine structure constant" a = e 2 /hvp ~ 1, 
since the Fermi velocity v F ~ 10 6 m/s w c/300 (c is the 
velocity of light), and therefore we have a strong-coupling 
version of QED. Below we set ft = c = 1. 

Existence of induced charge density in the electric field 
of heavy nucleus due to vacuum polarization is one of 
the most important effects of QED. This problem was 
inve stig ated in detail in many papers, see, e.g., Rcfs. 



13 . |20L l2ll , [22j . Charged impurity screening in graphene 



can also be treated in terms of vacuum polarization 
if 



25, m, 27, \m, \m \m 31, 32, & 



In Ref.& 



the induced charge density in graphene has been inves- 
tigated analytically using convenient integral represen- 
tation for the Green's function of the two-dimensional 
Dirac equation of electron in a Coulomb field. Calcula- 
tion of the induced charge has been performed exactly in 
the charge of impurity. In Ref. [32| , the Green's function 
has been obtained following the method based on the op- 



erator technique suggested in Ref. 34] . In the present pa- 
per, we use similar integral representation for the Green's 
function to derive the induced current in the field of in- 
finitesimally thin solenoid. Calculation is performed for 
arbitrary value of the magnetic flux 

The induced density and induced current in the vector 
potential 



A(r) 



$[i/xr] 
27rr 2 



where u is the unit vector directed along z-axis, have the 
form 

de. 



Pind 



(r) = -ieN 



2vr 
de. 



Tr{G(r,r|e)}, 



J ind (r) = -ieNvp — Tr{<xG(r, r|e)} , (1) 
Jc 2tt 

where N = 4 reflects the spin and valley degeneracies, 
and the Green's function G(r,r'\e) satisfies the equation 

[e - v F cr ■ (p - eA(r))] G(r, r'\e) = S(r - r')I. (2) 

Here <r = (cri,(72), and cti are the Pauli matrices; p = 
iPx,Py) is the momentum operator, r — (x,y), and 
/ = diagjl, 1}. The matrixes cr do not act on the spin 
variables but on the pseudo-spin ones and the spin de- 
grees of freedom are taken into account in a factor N. 
According to the Feynman rules, the contour of integra- 
tion over e goes below the real axis in the left half-plane 
and above the real axis in the right half-plane of the 
complex e plane. It is convenient to write the function 
G(r, r'\e) as 

G(r, r'|e) = [e + v F a ■ (p - eA)] D{r, r'\e) , (3) 

Where D(r,r'\e) is the Green's function of the squared 
Dirac equation, 

[e 2 - v%(p - eA(r)) 2 + v F $5(r)a 3 ] D(r, r'\e) 

= S(r - r')I. (4) 
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For r ^ and r' ^ , we can omit the term with 8- 
function so that 

[e 2 - vp(p - eA(r)) 2 ] D(r, r'\e) = S(r - r')I. (5) 

The equation §E§ has regular and singular solutions at 
r = and r' = 0. The Green's function of the Dirac 
equation with the magnetic-solenoid field was considered 
in Ref. [HI taking into account both regular and singular 
parts. The singular solutions originate from the singular 
behavior at r = of the vector potential A(r). There- 
fore, to find the correct superposition of regular and sin- 
gular solutions it is necessary to perform the appropriate 
regularization. If we take in mind a real solenoid, then 
the natural regularization is the finite radius R of this 
solenoid, see Ref.0]. Then it is possible to show that, to 
calculate the induced charge density and the induced cur- 
rent at r 3> R, we can use the regular Green's function 
of the equation §5$) while the singular part of the Green's 
function as well as the term with <5-function in Eq. (j4|) 
determine these quantities at r < R. The induced cur- 
rent and the induced charge density at r < R depend on 
the magnetic field distribution inside the solenoid and, 
therefore, are model-dependent. At r ^> R, the contri- 
bution of the singular part of the Green's function to the 
integrand in Eq.© contains a factor (eR)^ with some 
positive (3. The main contribution to the integral over e 
at r ^ R is given by the region e ~ 1/r so that the con- 
tribution of the singular part of the Green's function to 
the induced current is suppressed by the factor (R/r) 13 . 
This situation is completely similar to the problem of a 
finite nuclear size at the calculation of the vacuum po- 
larization effects in heavy atoms, see, e.g., Ref. [36| . In 
quantum field theory it is possible to consider rcgular- 
izations different from the solenoid radius R. As a result 
some uncertainty appears in the prediction of physical 
quantities 0, 0, Q ■ In the present paper we consider the 
induced charge and induced density at r ^ which are 
model- independent. Substituting the function D(r,r'\e) 
to Eq.© in the form 



D(r,r , |e) = — V e^-^ A m (r, r'\e)I , 

Z7T i — * 



(6) 
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and righthand side of the Eq.© as 

pf i\ oo 

o(r — r 



5(r~r')I 



2~K*Jrr 
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we obtain that the function A m (r, r'|e) satisfies the equa- 
tions: 



4 + ~4 - { -^) ^(r, r'\e) = %^ ,(7) 
vp r Or Or r z / v F yrr' 



where 7 = e$/(27r). One can see that the function 
A m {r, r'\e) can be obtained from the Green's function for 



the free radial Schrodingcr equation in the 2D case by the 
substitution m — ► m — 7. A convenient integral represen- 
tation for the function A m (r, r'\e) can be obtained using 
the operator method developed in Ref. [34[ at the calcu- 
lation of the Green's function for the Dirac equation of 
an electron in a Coulomb field in 3D space. This method 
was recently used in Ref. for the case of 2D space. 
It follows from the results of Ref. at zero Coulomb 
field that 



A m {r, r'\e) 



ds 



vp sinh s 



I 2E\/rr 1 \ 

x exp\iE(r + r') coths - iir\] J 2 \ : , (8) 

\ smh s / 



where E = c/vf, A = |m — 7I, J2\{x) is the Bessel func- 
tion, fj, = +1 if RcE > and n = -1 if ReE < . The 
sign fi takes into account the analytical properties of the 
Green's function. 

Taking into account the analytical properties of the 
Green's function, the contour of integration with respect 
to e can be deformed to coincide with the imaginary axis. 
After these transformations, we obtain that pi n d(r) = 
as a result of integration over e. This fact can be eas- 
ily explained because, due to the Furry theorem, pi n d(r) 
should be the odd function of 7 = e$/(27r). However, 
in this case we would obtain that pi„d(f) is pseudoscalar 
that contradicts to the parity conservation of the mass- 
less 2D Dirac equation. For Jmd(r) we have 



Jind(r) = -^§rW ^^^{m-^) 

m=0 



00 00 



f dE [ . d f exp[-2£rcoths]/ 2 A 
J J smh s 



x J dE 





/ 2Er 
\ sinh s 



,(9) 



where l2\(x) is the modified Bessel function of the first 
kind. We note that Ji n d( r )> Eq. ©, is an odd function 
of 7, in accordance with the Furry theorem. To have a 
possibility to change the order of summation and integra- 
tion, we introduce some quantity S <C 1 as a lower limit 
of integration over s. After that we take the integral over 
E and obtain 



2?r 2 r 



m— — 00 



00 

/ds . . 

—— — exp — 2As 
smh s 



(10) 



As should be, Jind(r) depends only on the fractional part 
7 of 7, I7I < 1. The quantity 7 is 7 = 7 — n for 7 > 
and 7 = 7 + n for 7 < 0, where n is a maximal integer 
number less than |7|. Then we perform summation over 
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m and set 5 = 0. Wc have 



_ , . cNvf 
Jvndir) = - n \v x r| 



27r 2 r 3 



<7.s 
sinh i 



7exp(-2|7|s) 



+7 exp(— s)- 



sh(27s) sinh(27s) 



sinh s 2 sinh 2 s 
Taking the integral over s we finally arrive at 



Jind{r) = 6 1 / F -F(7)curl (- 
107T V r 



(11) 



F(7) = (l-2|7|) 2 tan(7r 7 ) 



(12) 



It is interesting that Ji n d(r) equals to zero at I7I = 1/2. 
This may be explained as follows. Due to invariancc 
of Jind{f) under the substitution 7^7 — 1 we have 
^(ItI) = ~ 1)) an( i due to the Furry theorem it 

should be F(j) = — F(— 7). From these two relations we 
obtain that F(±l/2) = 0. 
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FIG. 1: Dependence of the function F(j), Eq. (Tl2)) . on the 
fractional part 7 of 7 = e$/(27r), where $ is the magnetic 
flux. 

The induced charge density and current in the presence 
of an infinitesimally thin solenoid were also considered in 
Refs. 



371 l38l |. The results of these papers contain the 



contribution of both singular and regular parts of the 
Green's function since a regularization for the field was 
not applied. Therefore, the expressions for the induced 
charge density and current contain uncertainty which 
does not allow one to make any explicit predictions for 
these quantities. Note that our results arc in agreement 
with the contribution of the regular part of the Green's 
function in Ref. [HI (first term in Eq.(6.14)). 

To summarize, we have investigated the effect of vac- 
uum polarization in the field of an infinitesimally thin 
solenoid at distances much larger than the radius of 
solenoid. It turns out that the induced charge density 
is zero. We have derived exactly in a magnetic flux the 
expression for the induced current. This current is a peri- 
odic function of the magnetic flux and is equal to zero not 
only at the integer values of 4>/<I>o but also at half-integer 
values of this ratio. Though the system considered in our 
paper consists of graphenc and a solenoid perpendicular 



to the plane of a sample, the results can be easily gener- 
alized for another systems such as studied in Ref. 9]. 
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